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Temporally-Correlated Channels 
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Abstract 

Feedback of quantized channel state information (CSI), called limited feedback, enables transmit 
beamforming in multiple-input-multiple-output (MIMO) wireless systems with a small amount of over- 
head. Due to its efficiency, beamforming with limited feedback has been adopted in several wireless 
communication standards. Prior work on limited feedback commonly adopts the block fading channel 
model where temporal correlation in wireless channels is neglected. This paper considers temporally- 
correlated channels and designs single-user transmit beamforming with limited feedback. Analytical 
results concerning CSI feedback are derived by modeling quantized CSI as a first-order finite-state 
Markov chain. These results include the source bit rate generated by time-varying quantized CSI, the 
required bit rate for a CSI feedback channel, and the effect of feedback delay. In particular, based on 
the theory of Markov chain convergence rate, feedback delay is proved to reduce the throughput gain 
due to CSI feedback at least exponentially. Furthermore, an algorithm is proposed for CSI feedback 
compression in time. Combining the results in this work leads to a new method for designing limited 
feedback beamforming as demonstrated by a design example. 

I. Introduction 

For a multiple-input-multiple-out (MIMO) communication system, transmit beamforming alleviates the 
negative effect of channel fading by exploiting spatial diversity, and thereby increases system throughput. 
Typically, transmit beamforming requires feedback of channel state information (CSI) from a receiver 
to a transmitter. Such CSI feedback can potentially incur excessive overhead due to the multiplicity 
of channel coefficients. This motivates designing highly-efficient CSI quantization algorithms based on 
communication measures such as information capacity. These algorithms for both beamforming and other 
MIMO techniques form an active field called limited feedback (see e.g. [1] and the references therein). 
This paper focuses on limited feedback beamforming systems over temporally correlated channels. It 
addresses a set of open issues including the information rate inherent in CSI, the required bit rate for 
a CSI feedback channel, the effect of feedback delay on throughput, and feedback CSI compression. 
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This paper provides a systematic method of designing the CSI feedback link of a limited feedback 
beamforming system. 

A. Prior Work and Motivation 

For simplicity, most prior work on limited feedback adopts the block fading channel model, where each 
channel realization remains constant in one block and different realizations are independent [2]. Using this 
model, designing limited feedback reduces to a vector quantization problem [3]. Different methods for 
quantizing CSI have been developed such as line packing [4], [5], combined parametrization and scalar 
quantization [6], subspace interpolation [7], and Lloyd's algorithm [8]. Furthermore, different types of 
limited feedback systems have been investigated, including beamforming [4], [5], precoded orthogonal 
space-time block codes [9], precoded spatial multiplexing [10], and multiuser downUnk [11]. Prior work 
focuses on analyzing and miiumizing CSI inaccuracy due to quantization. The existing results do not 
account for how channel coherence time influences the information rate generated by time-varying CSI, 
and the bit rate that a feedback channel must support. These issues are addressed in this paper 

For practical systems, the block fading assumption in prior work is pessimistic since channel tem- 
poral correlation often exists and can be used for compressing feedback CSI by incremental feedback. 
Algorithms for feedback compression have been proposed in [6], [12], [13] that exploit channel temporal 
correlation. In [6], a MIMO channel is parameterized and the feedback of each parameter is compressed to 
be one bit. Nevertheless, the multiplicity of the channel parameters compromises the feedback efficiency. 
In [12], the feedback CSI is compressed to be one bit but requires the transmitter to periodically broadcast 
channel subspace matrices. Building on the preliminary results of the current work in [14], variable-length 
source codes such as a Huffman code [15] are applied for CSI feedback compression in [13]. Despite 
minimizing the average CSI source bit rate, this approach may not suit practical applications where CSI 
feedback consists of fixed-length bit blocks [16], [17]. 

In a practical system with CSI feedback, feedback delay exists due to sources such as signal processing, 
propagation and channel-access protocols. The negative effects of CSI feedback delay on bit-error rate or 
information capacity have been observed in the literature of MIMO communication [18]-[24]. Specifically, 
this delay is found to decrease received signal power and cause interference between spatial data streams. 
The negative impact of CSI feedback delay can be alleviated to some extent by channel prediction [24], 
[25]. Despite its importance for designing MIMO systems, there exists no simple relationship between 
CSI feedback delay and throughput. This motivates the current work on deriving such a relationship in the 
context of limited feedback beamforming. The results of this work have been vaUdated by measurement 
data from a MIMO prototype over an indoor channel [26]. 
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The main approach of this work is to model quantized CSI as a finite-state Markov chain (FSMC), 
which allows the use of Markov chain theory as an analytical tool. A similar approach is common 
in modeling single-input-single-output (SISO) fading channels [27]-[31]. A FSMC model for fading 
captures the stochastic feature of wireless channels that is missing in block fading. Furthermore, the 
FSMC model is simple enough for allowing tractable performance analysis of communication systems. 
For these reasons, since it was proposed in [27], this model has be widely applied in wireless channel 
modeling, communication and networking. The FSMC models have been proposed for the satellite [32], 
indoor [33], Rayleigh fading [27], [29], [34], Nakagamni fading [35], and Rician fading [28] channels. The 
accuracy of FSMC models has been verified based on different criteria, including information capacity 
[36], [37], packet errors [38], burst errors [39], and autocorrelation functions [28], [29]. Due to its accuracy 
and simplicity, the first-order FSMC channel models have been used in designing and analyzing adaptive 
video encoding [40], maximum a posteriori decoding [41], downlink rate control [42], downhnk power 
control [43], automatic repeat request (ARQ) [44], and maximum likelihood detection [45]. FSMCs have 
been also adopted for modeling scalar functions of MIMO channel matrices such as the largest singular 
value [46], the condition number [47], and the Frobenius norm [48]. Nevertheless, there exist few results 
on using FSMC for designing MIMO limited feedback systems despite their popularity. 

B. Contributions 

This paper presents a set of analytical results useful for understanding and designing MIMO hmited 
feedback beamforming over a temporally-correlated channel. These results assume that the time-varying 
quantized CSI follows a first-order finite-state Markov chain, called channel state Markov chain. The 
contributions of this paper are summarized as follows. 

1) Quantized CSI is treated as a Markov data source and the corresponding information bit rate, called 
CSI source bit rate, is derived in terms of the probabilities of the channel state Markov chain. In 
particular, the CSI source bit rate is shown to decrease linearly with the probability that the channel 
state remains unchanged over two consecutive data symbol durations. Note that this probabihty is 
larger for longer channel coherence time and vise versa. The CSI source bit rate provides a measure 
of the rate of information generated by the temporally-correlated channel. 

2) CSI feedback can rely on periodic or random-access [49] feedback protocols, corresponding to 
separated or shared feedback channels for different users. Based on the periodic feedback protocol, 
the bit rate supported by a feedback channel, called CSI feedback bit rate, is derived under a 
feedback outage constraint, where an outage refers to multiple channel-state transitions within one 
feedback interval. This constraint guarantees that CSI feedback is sufficiently frequent. The derived 
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CSI feedback bit rate is useful for allocating bandwidth for the CSI feedback channel. 

3) Define the feedback throughput gain as the difference in throughput between the cases of delayed 
CSI feedback and no feedback. Based on the theory of Markov chain convergence rate, the feedback 
throughput gain is shown to decrease at least exponentially with feedback delay and inversely with 
the feedback interval. The exponential rate increases with decreasing channel coherence time and 
vise versa. This result enables the joint design of the mobility speed, the bandwidth of the feedback 
channel, and feedback delay sources such as signal processing complexity and the propagation 
distance, under a constraint on the feedback throughput gain. Moreover generally, the feedback 
throughput gains are proved to be delay sensitive. 

4) Finally, an algorithm is proposed for compressing feedback CSI by exploiting residual temporal 
correlation in feedback CSI. Note that the temporal correlation in CSI is largely reduced by the 
periodic feedback protocol. This algorithm compresses feedback CSI by truncating low-probability 
transitions between states of the channel state Markov chain. Moreover, this algorithm alternates 
compressed and uncompressed CSI feedback to prevent propagation of CSI errors due to such 
truncation. 

The differences between this manuscript and its conference versions [14], [50], [51] are highlighted as 
follows. First, the CSI feedback bit rate derived in this manuscript targets a periodic feedback protocol 
under a feedback outage constraint, but that in [14] is based on an aperiodic feedback protocol without 
such a constraint. Note that the periodic feedback protocol is more suitable for existing communication 
standards such as 3GPP-LTE [16], where a fixed number of bits in each reverse-link data block are 
allocated for CSI feedback. Second, this manuscript considers both fixed feedback delay and variable 
delay inherent in the periodic feedback protocol, but only the former is addressed in [50]. Moreover, the 
analysis on the effect of fixed feedback in [50] has been made more rigorous in this manuscript. Third, 
the accuracy of the Markov chain model for CSI feedback is verified by simulation in this manuscript 
but not in [14], [50], [51]. Fourth, the CSI feedback compression algorithm proposed in [51] is improved 
in this manuscript to avoid propagation of CSI errors due to truncation of low-probability channel-state 
transitions. Finally, a new design example is included in this manuscript for demonstrating the joint 
application of different analytical results from this work. 

C. Organization and flotation 

The remainder of this paper is organized as follows. In Section |II1 a limited feedback beamforming 
system is described. Section |lll] presents the definition and construction procedure of the channel state 
Markov chain. In Section |IVl CSI source and feedback bit rates are derived. Section |V] focuses on the 
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relationship between the feedback throughput gain and feedback delay. In Section |VT1 an algorithm for 
CSI feedback compression is proposed and analyzed, followed by concluding remarks in Section IVIIII 
Notation: Capitalized and small boldface letters denote matrices and vectors, respectively. The super- 
script * represents the complex conjugate and transpose matrix operation. The operator [•]£ gives the Ah 
component of a matrix. Similarly, returns the m)th component of a vector. 

II. System Description 

The limited feedback beamforming system illustrated in Fig. [U can be separated into the forward and 
the CSI feedback links, which are described in the following sub-sections. In this system, all signals are 
discrete-time and sampled at the sampling rate 1/T, where T denotes one sample interval. Without loss 
of generality, T is set equal to one symbol duration. The sample index is denoted by the subscript n. 

A. Forward Link 

The forward link refers to the data path in Fig. [T] from the input of the beamformer to the output of 
the maximum-ratio combiner. The nth received data symbol after maximum ratio combining, denoted as 
y„, is given as 

Vn = \/Pw*H„f„a;„ + i/„ (1) 

where x„ and i/„ are CM{0, 1) random variables modeling respectively the nth transmitted data symbol 
and the nth sample of the additive white Gaussian noise (AWGN) process, and P is the transmit signal- 
to-noise ratio (SNR). Let Nt and Nr denote the number of transmit and receive antennas, respectively. 
Then the Nr x 1 complex unitary vector w„ represents the weights for maximum ratio combining, and 
the Nt X 1 complex unitary vector f„ denotes the transmit beamforming vector. The Nr x Nt matrix 
represents the nth realization of the MIMO channel. The random process {H„} is assumed to be 
stationary and temporally correlated. Note that the common assumption of a complex Gaussian channel 
is unnecessary for our analysis. 

The receiver continuously estimates the CSI sequence {H„} using pilot symbols sent by the transmitter. 
The estimated CSI is used for computing the maximum-ratio combining vector and the beamforming 
vector for feedback. This paper considers the scenario where CSI quantization and feedback delay are 
the main sources of transmit CSI inaccuracy. Thus, CSI estimation is assumed perfect for simplicity. 
This assumption is commonly made in the literature of limited feedback (see e.g. [4], [5]). Based on 
this assumption, the maximum-ratio combining vector is computed as w„ = H„f„/||H„f„||. In the next 
section, the beamforming vector f„ is selected from a codebook for maximizing the receive SNR. 
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Fig. 1. Limited feedback beamforming system 
B. CSI Feedback Link 

In Fig. [U CSI feedback link refers to the CSI path from the input of CSI quantizer to the output of 
the feedback channel. To satisfy a finite -rate feedback constraint, CSI is quantized efficiently by using a 
Grassmannian codebook, denoted as J^, designed for beamforming in [4], [5]. Let denote the ith of 
N unitary vectors in J^. To maximize the receive SNR, the quantizer function, denoted as (or Qj), 
maps the channel matrix H„ to a beamforming vector in the codebook (or its index denoted as In) 
as follows [4], [5] 



Q/(H„) = argmax||H„v||^ , /„ = Qi(H^ 



arg max llHnV/'l 



(2) 



The feedback of the index /„, called the channel state hereafter, is sufficient for the transmitter to retrieve 
the selected beamforming vector f„ from the codebook J^. 

From Q, the channel states, aie a sequence of alphabets of N letters. This alphabet sequence 

is encoded by using a B-hit fixed-length code, where B = \og2 N [15]. The CSI bit rate at the output of 
the source encoder is derived in Section [TV-AI Intuitively, a more efficient alternative is to encode a long 
block of channel states by using variable-length codes such as a Huffman code [15]. Nevertheless, block 
CSI encoding contributes additional feedback delay, which decreases throughput significantly as shown 
in Section jV] Furthermore, a variable CSI codeword length is unsuitable for typical limited feedback 
systems such as 3GPP-LTE [16], where an uniform number of bits are allocated for each instant of 
periodic CSI feedback. This also motivates the use of a periodic feedback protocol in this paper. 

The periodic feedback protocol transmits the latest encoded CSI to the feedback channel at fixed 
intervals. Thereby this protocol introduces variable CSI feedback delay as elaborated shortly. Let the 
integers K and D denote respectively a feedback interval and fixed feedback delay in samples, where D 
is contributed by sources including signal processing and propagation. Then CSI used for beamforming in 
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successive symbol durations lags behind the corresponding channel states by D, (D + 1), • • • , {D^K—X) 
samples repeated in a cyclic order. The CSI feedback bit rate for the periodic feedback protocol is derived 
under a feedback outage constraint in Section [TV-Bi and the effects of CSI feedback delay are investigated 
in Section IVl 

To simplify analysis, the feedback channel is assumed free of errors, which is typical in the literature of 
limited feedback (see e.g. [1], [4], [5]). This assumption is justified by the fact that beamforming feedback 
as a control signal is usually well protected by using error-correction coding or high transmission power. 

III. Channel State Markov Chain 

In this section, the channel state sequence is modeled as a first-order finite-state Markov chain. 
The channel state Markov chain is used as the primary analytical tool in the sequel. As discussed in 
Section Jl finite-state Markov chain models for fading have been established as valid models of wireless 
channels [27]-[32], [34]-[39], and widely used in wireless communication and networking [40]-[45]. 
Following this common approach, a finite-state Markov chain is used for modeling partial CSI of a 
MIMO fading channel in this section. The resultant Markov chain is called the channel state Markov 
chain. Its parameters, namely the state space and the stationary and transition probabilities, are related to 
the channel statistics and the CSI quantizer codebook. The channel state Markov chain model is validated 
using simulation in Section IVII-AI 

The channel state Markov chain models the time-variation of the channel state I„ in Q. Mapped from 
a stationary channel by the quantizer function in Q, {/„} is a finite-sate stationary stochastic process. 
Motivated by the common approach in fading channel modeling [27]-[31], this process is modeled as a 
homogeneous finite-state Markov chain of order one with the state space T = {1, 2, • • • , A^} [52]. The 
key property of this Markov model is that conditioned on the most recent state ln-\. In is independent 
of the past states In-2, In~3, ■ ■ ■ , thus 

Pr(/„ = in I In-1 = tn-lJn-2 = 4-2, • • • ) = = 4 | In-l = In-x)- 0) 

Give the above property, the probability of a transition from state m to £ is defined as Pim = PrC-^n = 
i I ln-\ = m). Moreover, the stationary probability of state I is defined as vr^ = Pr(/„ = £). For 
convenience, denote the x 1 stationary probability vector as tt with = vr^, and the N x N 
transition probability matrix as P with [P]em = Pem, which is known as the stochastic matrix [52]. 

The channel state Markov chain is assumed to be ergodic [52]. The ergodicity implies three properties. 
First, the states of the Markov chain are communicating, namely that a transition between every pair of 
states occurs within a finite duration. Second, each state is recurrent and thus the probability of returning 
to a same state is one. Third, each state is aperiodic. This property exists if all possible time durations 
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(in samples) of leaving and returning to each state have the common divisor of one. The assumption of 
ergodicity for the channel state Markov chain is valid for typical continuous channel distributions such 
as Rayleigh or Rician [2]. This assumption provides the following property [52] 

lim = [7r,7r,--- ,7r]. (4) 

The probabilities vr^ and Pg^ depend on both the channel statistics and the CSI quantizer codebook. 
To specify these relationships, define a Voronoi cell based on the quantizer in ([2]) as [3] 

= {H e C^'-x^' I llHv^ll > ||Hv„|| V m G X, m / ^} (5) 

where denotes the £th unitary vector in the codebook J^. This set Ve maps the channel to the ^th state 
of the channel state Markov chain as follows 

I^=£^Hne Ve. (6) 

Using the above relationship, the probabilities of the channel state Markov chain can be related to the 
channel statistics as 

TTe = Pr (H e Ve) , Pem = Pr (H„ G Ve \ H„_i G V„) . (7) 

In general, dV]) does not yield closed-form expression for the Markov chain probabilities except for the 
degenerate case of single antennas. Nevertheless, these equations are useful for computing the probabilities 
by simulation as discussed in the next section. 

In this paper, the channels and hence the channel state Markov chains are assumed stationary. It is 
assumed that the receiver perfectly estimates the Markov-chain parameters, namely the stationary and 
transition probabilities. To facilitate CSI feedback, the receiver communicates once to the transmitter 
functions of the Markov chain parameters including the CSI feedback bit rate, the allowable feedback 
delay, and the source codes for compressing feedback CSI, which are derived in the sequel. Given channels 
that are stationary or approximately so, the infrequent feedback of these functions incurs overhead that 
is much smaller than that for the instantaneous CSI feedback for beamforming, and hence is neglected 
in our analysis. 

IV. CSI Source AND Feedback Bit Rates 

In this section, the overhead required for CSI feedback is analyzed based on the channel state Markov 
chain. Specifically, the CSI source and feedback bit rates are derived in Section IIV-AI and Section IIV-BI 
respectively. 
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A. CSl Source Bit Rate 

The CSI source bit rate, defined as tiie average rate of encoded CSI bits, measures the amount of 
information generated by the temporally correlated MIMO channel. Based on fixed-length source coding 
as discussed in Section |II1 the CSI source rate is the product of the CSI codeword length and the average 
rate of channel-state transition given as 



Given the ergodicity assumption for the channel state Markov chain, Rg in ([D can be rewritten in terms 
of the Markov chain probabilities by applying the ergodic theorem in Lemma [T] from [53, Corolary 4.1]. 

Lemma 1: Given a function z : I x I ^ M, that satisfies X]m=i k(^) "^)|^^-P£m < oo, the 

following convergence exists almost surely 



The main result of this section as summarized in Proposition [Tf follows from ([8]) and Lemma [T] Note that 
the condition for applying Lemma [T] is checked to be satisfied by setting z{i,m) = 1{1 / m). 
Proposition 1: The CSI source bit rate is 



The above result shows that Rg decreases linearly with increasing probabilities {Pu}, which characterize 
the degree of channel temporal correlation. By definition, P££ is the probability that the channel state 
remains as £ for two consecutive samples. The values of {Pee} are close to one if the channel is highly 
correlated in time, or close to zero for fast fading. As shown by simulation in Section [VII-Bi Rg increases 
linearly with Doppler shift for spatially i.i.d. Rayleigh fading and the Clarke's fading model [2]. 

B. CSI Feedback Bit Rate 

The CSI feedback bit rate is defined as the maximum bit rate supported by the feedback channel. 
In this section, the CSI feedback bit rate is derived under a constraint on feedback outage probability. 
Recall that a feedback outage refers to the event that more than one channel-state transition occurs within 
one feedback interval. Let Z denote the number of channel-state transitions in a feedback interval of K 
samples. Moreover, let Pq represent the feedback outage probability, thus Pq := Pr(Z > 1). For a small 
real number 6 > 0, the constraint Pq < S ensures that CSI feedback is sufficiently frequent. 

The main result of this section is obtained and summarized in the in following proposition. 




(8) 



n=l 




(9) 




(10) 
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Proposition 2: Under the feedback outage constraint Pq < 5, the CSI feedback bit rate is Rj = 
bit/s where the feedback interval in samples, K, is the largest integer subject to 

N N N pK^p ^ 

E ^^^"^ + E E p 1 p ^ 1 - ^- (11) 

1=1 1=1 m=l "^"^ 

Proof: See Appendix lAl □ 

To observe the dependence of the CSI feedback bit rate on channel temporal correlation, bounds on Rj 
can be obtained as a by-product of the proof for Corollary [T] 

log(maXa Paa) ^ r> ^ 

log(l -6)- log (l + Em^, p,^) - l°g(l-^) 

The numerators of the above bounds suggest that Rf decreases with increasing probabilities {Pee}, which 
characterize the degree of channel temporal correlation as mentioned earlier. In other words, Rj is smaller 
for a more temporally correlated channel and vise versa, which is also observed for the CSI source bit 
rate given in Proposition [T] 

The feedback interval K for the extreme cases of small and large target feedback outage probabilities 
are characterized in the following corollary of Proposition |2l 

Corollary 1: For a small target feedback outage probability 5 — > 0, the normalized feedback interval 
converges to one: K 1; for a large probability 5 ^ 1, K scales as 

K = ri\og{l-5) + 0{l) (13) 

where ^ — 5- < n < 1 — s-. 

log maxf Fte — ' — log miiif Fie 

Proof: See Appendix iBl □ 

For (5^1, the first term ?7log(l — 5) in (fT3] ) is dominant because it is asymptotically large. The factor 
r] in ( fT3] ) represents the sensitivity of K to the change on the feedback constraint 5. This factor depends 
on the degree of channel temporal correlation through the probabilities {Pu}- 

V. Feedback Delay 

This section focuses on the effects of feedback delay on the throughput of the limited feedback 
beamforming system. We consider hoih fixed delay due to signal processing and propagation, and variable 
delay caused by the periodic feedback protocol. In Section IV-A[ the feedback throughput gain with fixed 
feedback delay is defined and derived. In Section |V-B[ fixed feedback delay is shown to reduce the 
feedback throughput gain at least at an exponential rate. Building on this result, an upper-bound on the 
feedback throughput gain is derived as a function of both fixed and variable feedback delay in Section lV-Cl 
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A. Feedback Throughput Gain 

The feedback throughput gain is defined as the gain in ergodic throughput due to delayed CSI feedback 
with respect to the case of infinite feedback delay. In this section, only fixed feedback delay is considered 
and a corresponding upper bound on the feedback throughput gain is derived. This result is extended 
in Section IV-CI to include the effect of variable feedback delay. Let D denote fixed feedback delay 
in samples, and R represent ergodic throughput. Then the feedback throughput gain can be written as 
Ai?(L>) = R(p) -R{oo). 

The ergodic throughput R{D) is derived as follows. With feedback delay, the input-output relationship 
for the limited feedback beamforming system in ([Hi can be rewritten as 



where the effective channel gain gn{D) = ||H„f and f{In-D) is the transmit beamforming 
vector chosen based on the delayed feedback CSI In-D- With only beamforming feedback, gn{D) is 
unknown to the transmitter, and hence constant transmission power P is optimal [2]. Given feedback 
delay, deriving the optimal strategy for transmit beamforming is difficult as it depends on the codebook, 
the channel stochastic distribution, and the channel prediction algorithm. For simplicity, this paper 
adopts the strategy of applying the codebook vector corresponding to In-D as the beamforming vector. 
Equivalently, i{In-D) = Q/(H„_/)) with the quantizer function Q{-) given in Q. Note that this sub- 
optimal beamforming strategy is optimal for zero feedback delay [4]. Using this strategy, the ergodic 
throughput of the effective scalar channel in (fT4l) is obtained as [2] 



To achieve this throughput, the channel code for forward-link data is assumed to be sufficiently long to 
attain channel ergodicity [2]. In other words, the code length and hence the decoding delay are much 
longer than channel coherence time. If a constraint on the decoding delay is required, the outage capacity 
is a more appropriate performance metric [54]. 

The ergodic throughput in ([T5] ) can be rewritten in terms of the probabilities of the channel state 
Markov chain as shown in the following lemma. 

Lemma 2: For fixed feedback delay of D symbol durations, the ergodic throughput is 



Un = \fPgn{D)Xn + f, 



■n 



(14) 



R{D) = E [log2(l + gl{D))] = E [log2(l + P||H„Qy(H„_B)f )] . 



(15) 



N 



N 




Im "™ 



(16) 



£=1 m=l 



where 



i?^^ = E[log2(l + P||H„v, 



m 



(17) 



Proof: See Appendix |C] 



□ 
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In general, the constants {Rim\ do not have closed-form expressions and have to be estimated by using 
Monte Carlo simulation. 

Consider the extreme cases of zero delay {D = 0) and infinite delay {D oo). The corresponding 
ergodic throughput is obtained and shown in the following corollary of Lemma |2j 

Corollary 2: The ergodic throughput for zero feedback delay R{0) and infinite delay R{oo) are 

N 

R{0) = J^i2«7r, = E[log2(l + P||HQ^(H)||2)] (18) 

N N N 

5^E[log2(l + P||Hv,||2)]7r,. (19) 

e=l m=l 1=1 

Proof: See Appendix IDI □ 
A few remarks are in order. 

1) For zero feedback delay, -R(O) is equal to the ergodic capacity for limited feedback beamforming 
based on the block fading channel model [4, (26)], where feedback delay is omitted. 

2) For infinite feedback delay, the beamforming vector at the transmitter is adapted to obsolete CSI 
independent of the current channel state. Higher ergodic throughput than R{oo) may be achieved 
by using space time block coding or adapting the beamforming vector to channel statistics [55], 
both require no instantaneous CSI. Despite its sub-optimality, R{oo) serves as a reference value 
for computing the feedback throughput gain in this paper. 

3) It can be observed from (fTSl) and ( fT9l ) that both i2(0) and R{oo) are independent of the transition 
probabilities of the channel state Markov channel. This suggests that channel dynamics affect only 
the case of finite feedback delay as considered in the next section. 

Finally, from Lemma [H Corollary |2] and the definition, the feedback throughput gain is readily written 
as shown in the following proposition. 

Proposition 3: Assuming only fixed feedback delay, the feedback throughput gain is given as 

N N 

Ai?(I)) = ^^i?,^{[P^]^^-7rj7r™ (20) 

e=i m=l 

where Rim i^ given in (117b . 

An upper-bound on R{D) is derived in the next section. 

Other metrics: Besides the feedback throughput gain, other metrics can be defined for quantifying the 
effects of CSI quantization and feedback delay on the throughput of a transmit beamforming system. The 
ergodic capacity for perfect CSI feedback is Cideai = IE [log2 (l + PAf (H))] where Ai(H) represents the 
largest singular value of H. With respect to the case of perfect CSI feedback, the capacity loss due to CSI 
quantization, called quantization loss, can be defined as VgC = Cideai — ^(0) where R{0) follows from 
([T5] ). As shown in [4], [5], quantization loss is small given just a few bits of resolution for quantized 



Huang et ah: Limited feedback beamforming over temporally-correlated channels 



13 



CSI. Moreover, the capacity loss due to both CSI quantization and feedback delay can be written as 
= Cideai — R{D). Next, using the case of infinite feedback delay as the reference, the maximum 
throughput gain is defined as Amax^ = Cideai — R{oo). The feedback throughput gain defined earlier as 

A. R{D) = R{D) — R{oo) takes into account of both CSI quantization and feedback delay. By definition, 
the feedback throughput gain IS.R{D) is the complement of the capacity loss VqdR in the sense that 

Ai?(D) + \dR = AmaxC. 

B. Effect of Fixed Feedback Delay 

In this section, increasing fixed feedback delay is shown to reduce the feedback throughput gain at 
least at an exponential rate. This result is derived based on theory of Markov chain convergence rate. To 
state this theory, define the time reversal of the stochastic matrix P, denoted by P, as a matrix whose 
components are given as [52] 

Pme = —Pim, l<m,i<N. (21) 

As implied by its name, P represents channel-state transitions in the reverse direction of P. For the 
special case of P = P, the channel state Markov chain is known as reversible [52]. Using the definition 
in ((2TI) . [56, Theorem 2.1] is restated as the following lemma. 

Lemma 3: For the ergodic channel state Markov chain, the following inequality holds 

lL\[^'']^m-M] <— ' ^<m<N (22) 

\£=1 / ™ 

where A G [0, 1] is the second largest eigenvalue of the matrix PP. 

By using Lemma [3l the main result of this section is derived and shown in the following proposition. 
Proposition 4: For feedback delay of D samples, the feedback throughput gain is upper bounded as 

< Ai?(D) < a(VA)-° (23) 

where a = X]m=i \/7rm max^ R^^ with R^^ given in ( fTTl ). and X G [0, 1] is the second largest eigenvalue 
of the matrix PP. 

Proof: See Appendix 10 □ 
A few further remarks are in order. 

1) The eigenvalue A is a key parameter characterizing the channel dynamics. A larger value of A 
indicates longer channel coherence time and vise versa. 

2) As observed from (1231 ). the feedback throughput gain decreases at least exponentially with the 
feedback delay. The decreasing rate is \/A and thus depends on channel coherence time. 
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3) For a reversible channel state Markov chain with P = P, \f\ is the second largest eigenvalue of 
the transition matrix P. 

4) As observed from simulation results in Section IVII-CI the upper bound in (1231) is tight in several 
cases of interest. 

C. Effects of both Fixed and Variable Feedback Delay 

In this section, the effects of both fixed and variable feedback delay on the feedback throughput gain 
are jointly considered. As discussed in Section III-BI due to the periodic feedback protocol, CSI used for 
beamforming in successive symbol durations encounters variable feedback delay of D, (D + 1), • • • , (-D + 
K — \) samples repeated in a cyclic order. Let R(K^ D) denote the ergodic throughput with both fixed 
and variable feedback delay. Thus, R{K, D) can be written as 

1 

R{K,D) = -Y,R{k + D) (24) 

fc=0 

where R{-) is given in (fT6l ). Define the corresponding feedback throughput gain as IS.R = R{K,D) — 
R{oo). Then we have the following corollary of Proposition |4l 

Corollary 3: For fixed feedback delay of D samples and a feedback interval of K samples, the 
feedback throughput gain is upper bounded as 

0< AR{K,D) < a{Vxf • (25) 

where a and A are identical to those in Proposition E] 

The upper-bound in (1251 ) can be decomposed into three factors a, (\/A)^ and ^(^) _ They characterize 

i\ (1— V A) 

the effects of CSI quantization, fixed feedback delay and the periodic feedback protocol, respectively. 
In particular, the feedback throughput gain decreases at least exponentially with D and approximately 
inversely with K. Thus, for small values of A corresponding to fast fading, the fixed feedback delay has 
a more significant effect on throughput than the delay due to the feedback protocol. 

Define the normalized feedback throughput gain as AR{K,D) = AR{K, D) / AR{0), where the 
normalization factor AR{0) corresponds to the ideal case of zero feedback delay. Motivated by the 
result in Corollary |3l AR{K, D) can be approximated as 

"l - (^/A)^ 



AR{K,D) w (VA) 



K{1 - VX) 

Id V 



(26) 



where the factors qu and qx separate the effects of the fixed and variable feedback delay on the feedback 
throughput gain. For a sanity check, the above approximated expression gives 1 for zero feedback delay 
K = D = and zero for infinite delay K ^ oo ov D ^ oo. The approximation of the normalized 
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feedback throughput gain in (l26l ) is observed to be accurate from simulation results in Section IVIII The 
result in (l26l) is useful for computing the allowable feedback delay under a constraint on the normalized 
feedback throughput gain as demonstrated by the design example in Section IVII-EI 



An algorithm is proposed in Section IVI-AI for compressing CSI feedback by exploiting temporal 
correlation in feedback CSI. The effect of feedback compression on feedback throughput gain is analyzed 
in Section IVllB] 

A. Algorithm 

The feedback compression algorithm is motivated by the existence of residual temporal-correlation 
in feedback CSI, which is largely reduced by the periodic feedback protocol (cf. Section IIV-BI ). To 
characterize such residual redundancy, feedback CSI is modeled as a Markov chain obtained by down- 
sampling the channel stat Markov chain at feedback intervals. Consequently, the stochastic matrix for 
this down-sampled Markov chain is with the {£,m)th component denoted as P^^^ ■ The correlation 
between two consecutive instants of feedback CSI is reflected in an uneven distribution of the transition 
probabilities in P^. Specifically, conditioned on the previous instant of feedback CSI, the current one 
belongs to a subset of the Markov chain state space with high probability. To define this subset, let the 
transition probabilities for the Markov chain state m, namely |-P£^^| , be indexed according to the 

(K) (K) (K) 

descending order of their values, thus -P^^ > " " " ^ ^Nm- ^'^^ ^ small positive real number e, the 
high-probability subset, called the e-neighborhood, of the mth Markov chain state is defined as 



This e-neighborhood groups most likely transitions from the Markov chain state m with total probability 
larger than (1 — e). High channel temporal correlation results in small e-neighborhoods and vise versa. 

Next, the feedback compression algorithm compresses feedback CSI in alternate feedback intervals to 
prevent propagation of CSI errors due to lossy compression. The possibility of error propagation is due 
to backward dependence of CSI compression in time. Conditioned on the nth feedback channel state 
IhK = that is uncompressed, the next one I(^n+i)K is compressed by lossy source coding. Specifically, 
I(n+i)K is encoded into one of |A/'m(e)| fixed-length codewords if I(n+i)K belongs to the e-neighborhood 
of InK, namely Mmif) defined in (|27] ). Otherwise, a codeword indicating CSI truncation is generated by 
the source encoder. It follows that compressing I(n+i)K requires a fixed-length source codebook having 
a size of [|A/'m(e)| + 1] and a codeword length of Bm(e) = riog2(|A/'m(e)| + 1)] bits. The above source 
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(27) 
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coding algorithm compresses CSI since Bm{e) is usually much smaller than B = log2 for the case 
of no compression. The number of bits for compressed CSI feedback should be uniform and thus is 
chosen as B{e) = max^i?£(e). Finally, to decode both compressed and uncompressed CSI feedback, the 
transmitter stores all (A^ + 1) source codebooks including one for uncompressed and for compressed 
CSI feedback. On decoding the codeword indicating CSI truncation, transmitter randomly chooses a 
beamforming vector from the codebook J^. 
A few remarks are in order. 

1) The proposed CSI compression algorithm alternates B and B{e) feedback bits for successive 
feedback instants, and hence the CSI compression ratio is 

A/i, = ^iM. (28) 

^ 2B 

This ratio is evaluated by simulation in Section IVII-DI 

2) The proposed algorithm can be also applied for compressing other types of feedback CSI such 
as channel-gain feedback useful for power control and scheduling. Different types of compressed 
feedback CSI can be integrated such that the total number of CSI bits per feedback instant remains 
constant, which suits practical systems such as 3GPP-LTE [16]. 

3) As mentioned earlier, the proposed feedback compression algorithm is preferred to the conventional 
lossless data compression algorithms, which use variable-length source coding and optionally 
block processing for higher compression efficiency [13], [15]. The reason is that variable-length 
CSI feedback is unsuitable for practical systems. Furthermore, block processing causes additional 
feedback delay that significantly decreases throughput as shown in Section IV-BI 

4) The proposed feedback compression algorithm complements the periodic feedback protocol on 
further decreasing the CSI feedback rate by exploiting channel temporal correlation. In contrast, 
based on an aperiodic feedback protocol, the original version of this algorithm presented in [51] is 
the only feedback compression function in the CSI feedback link. 

5) In practice, the receiver estimates the parameters of the channel-state Markov chain and computes 
the source codebooks for compressing feedback CSI. These codebooks are sent to the transmitter 
for decoding the compressed CSI. For stationary channels, one-time codebook feedback is sufficient. 

B. Feedback Throughput Gain with Feedback Compression 

The CSI feedback compression algorithm proposed in the preceding section is lossy and hence in- 
evitably incurs loss on the feedback throughput gain. To characterize this loss, the ergodic throughput 
with feedback compression, denoted as R{K,D), is written as 

R{K, D) = ^R{K, D) + ^1^R[K, D) + |fi(oo) (29) 
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where R(K^ D) is the ergodic throughput without feedback compression. The three terms on the right- 
hand-side of ( [29l ) correspond to uncompressed, compressed, and truncated feedback CSI, respectively. 
Note that the probability for CSI truncation is e. The feedback throughput gain with feedback compression 
is obtained by using ( |29l ) and given in the following proposition. Also shown is the CSI feedback bit 
rate for CSI compression. 

Proposition 5: The proposed CSI feedback compression algorithm has the following properties. 

1) The feedback throughput gain is 

AR{K, D, e) = (l - I) AR{K, D) (30) 

where AR{K, D) corresponds to the case of no feedback compression. 

2) The average CSI feedback bit rate is 

^ , , B + B(e) 

fl/W = (31) 
From (l30l ). e should be small so as to minimize the throughput loss due to lossy feedback CSI compression. 
Nevertheless, too small e decreases compression efficiency since -R/(e) converges to the feedback bit rate 
of uncompressed CSI, Rf, as e reduces to zero. According to (OTT ). Rf{e) is at least half of Rf. 

C. Extension: Block Feedback Compression 

For the feedback compression algorithm proposed in Section |VI-A[ compressing every other feedback 
CSI instant is usually sufficient since temporal correlation in CSI is substantially reduced by the periodic 
feedback protocol. Nevertheless, if the feedback interval is small, the residual correlation is significant 
and can be fully exploited by compressing CSI for multiple consecutive feedback instants, called block 
feedback compression. The algorithm in Section IVI-AI can be extended to block feedback compression 
as follows. Let W denote the compression block length in feedback interval KT. Feedback compression 
specified in Section IVI-AI repeats for a block of W feedback instants. Upon the first occurrence of 
feedback truncation in compressing a CSI block, the codewords indicating truncation are sent for the 
remaining feedback instants in this block. Uncompressed CSI feedback is performed between every two 
blocks of compressed CSI feedback to stop propagation of CSI errors due to feedback truncation. 

For block feedback compression, the ergodic throughput in ( |29l ) is rewritten as 
1 I ^ 

^) = wTT^^^' + wTi ^ { " ^) + [i - (1 - } • (32) 

Using ( [32l) . the feedback throughput gain is obtained as shown in the following lemma. 

Lemma 4: For block feedback compression, the feedback throughput gain is 

1 _ (1 _ 

Ai?(K, Z), e) = , Ai?(K, D) (33) 
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and the average CSI feedback bit rate is 

^f^'^ - {E + l)KT- ^ ^^^^ 

As observed from (l33l) . if the temporal correlation in feedback CSI is high (e << 1), AR{K, D, e) 

AR{K, D) and hence block feedback compression has a negligible effect on feedback throughput gain. 

Otherwise, the block length should be set as = 1 to avoid throughput loss. According to ( [34l ). the 

block feedback compression ratio can be up to \/{E + 1) instead of 50% for the original algorithm in 

Section lyTAl 

VII. Numerical Results and Design Example 

In this section, numerical results and a design example are provided to obtain further insights into 
designing limited feedback beamforming systems. The channel state Markov chain is constructed based 
on the spatially i.i.d. Rayleigh fading 4x4 MIMO channel model, corresponding to rich scattering and 
the size of antenna array equal to four. In this model, each channel coefficient is CAA(0, 1). The temporal 
correlation of each channel coefficient is assumed to follow Clarke's model [2] and is characterized by 
the zeroth order Bessel function of the first kind J'Q{2irf£)T), where //) is Doppler shift and r is the 
time difference between two samples of a channel coefficient [2]. Based on the above channel model, a 
first-order channel state Markov chain is constructed by simulation. 

A. Accuracy of Channel State Markov Chain 

In this section, the autocorrelation of quantized CSI is used as the criterion for validating the first-order 
channel state Markov chain. It is worth mentioning that the accuracy of this model is also supported 
by measurement data [26], where the feedback throughput gain due to transmit beamforming has been 
measured by using a MIMO-OFDM prototype over an indoor wireless channel. The measurement data 
closely matches the theoretical result in ( |26l ) derived using the channel state Markov chain. 

The quantized CSI for transmit beamforming is a correlated sequence of unitary vectors {fn}- From the 
quantizer functions in ([2]), it is clear that the CSI sequence f„ is one-to-one mapped to the channel state se- 



lf*f -I 

\'-n'-n+T\ 



where 



quence {In}- The autocorrelation function of the CSI sequence is defined as A{f) = E 
f is time separation in sample. The autocorrelation function can be evaluated by 1) simulating the 4x4 
MIMO channel as described earlier, 2) generating the quantized CSI sequence by applying dJ]), and 3) com- 
puting autocorrelation function based on channel ergodicity, namely A{f) = lim Yl^=i \^n^n+i-\^ ■ 
The autocorrelation function such evaluated is plotted in Fig. [2] with the legend "Clarke's model". The 
normalized Doppler shift is f^T = 10^^. By proper scaling the time separation (x-axis). Fig. |2]is also 
applicable for other normalized Doppler shifts. 
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Fig. 2. CSI autocorrelation versus time separation for the normalized Doppler shift /oT — 10~^. For the "Clarke's model", the 
autocorrelation function of each channel coefficient is given by the Clarke's function; for the "Markov chain", the autocorrelation 
of quantized CSI is modeled by a first-order Markov chain. 

Based on the channel state Markov chain, the CSI autocorrelation function is given by 



where is the fth vector in the codebook T . The above autocorrelation is plotted in Fig. |2] and labeled 
by the legend "Markov chain". 

As observed from Fig. |2j the first-order Markov chain provides an good approximation of the Clarke's 
fading model. In particular, the approximation is very accurate for correlation larger than 0.5. Moreover, 
both curves in Fig. [2] converge as the time separation increases. Fig. [2] suggests that higher-order channel- 
state Markov chains practically provide no improvement on modeling accuracy. For instance, two curves 
deviate significantly only for time separation larger than 10 samples. Consequently, a Markov chain of 
an order larger than 10 is required to provide an significantly better approximation of Clarke's model. 
Unfortunately, such a high-order Markov chain is impractical as its complexity increases exponentially 
with the order. 

B. CSI Source and Feedback Bit Rates 

Fig. |3] plots the normalized CSI source bit rate R^T in (fTOl ) versus normalized Doppler shift f^T on 
a logarithmic scale. The size of the CSI codebook is = {64, 128, 256}. As observed from Fig. |3j RsT 
increases Unearly with /dT. Moreover, the increasing rate differs for different values of A^ as reflected 
on the relative shifts between curves in Fig. |3] This rate is higher for larger A^ and vise versa. These 
agree with the intuition that faster fading and finer CSI quantization generate CSI at a higher rate. 
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10"* 10"' 10"' 

Normalized Doppler frequency 

Fig. 3. Normalized CSI source bit rate versus normalized Doppler shift for the CSI codebook size A'^ = {64, 128, 256}. 




Normalized feedback bandwidtfi 

Fig. 4. Feedback outage probability versus the CSI feedback bit rate normalized by the CSI source bit rate; the normalized 
Doppler shift is /dT = {10"*, 2 x 10"*, 3 x 10"*, 10"^}, and the CSI codebook size is iV = 128. 

Fig. |4] plots feedback outage probability Po versus the normalized CSI feedback bit rate Rf / R^. The 
normalized Doppler shift is foT = {10"'*, 2 x 10-^,3 x 10-^,10-3}, and the CSI codebook size is 
N = 128. As observed from Fig. |4l for Rf/Rs of practical interest (e.g. smaller than 6), Pq is insensitive 
to a change on foT. As observed from Fig. |4l Rf/Rs must be larger than 2 to keep Po reasonably low 
(e.g. smaller than 0.1). By using Fig. |4] and Fig. [3l the CSI feedback bit rate can be obtained for a target 
feedback outage probability and a Doppler shift. 
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Fig. 5. Comparison of the effects of CSI quantization and feedback delay on the ergodic throughput of a transmit-beamforming 
system. Three cases of CSI feedback are considered, namely 1) perfect, 2) quantized, and 3) quantized and delayed CSI feedback. 
The number of antennas in an array is L = 4, the SNR is 10 dB, the normalized Doppler shift is /uT = 10~^, and the quantizer 
codebook size is TV = 128. 

C. Feedback Delay and Feedback Throughput Gain 

Fig. [5] compares the effects of CSI quantization and feedback delay on the ergodic throughput of a 
beamforming system. To be specific, the ergodic throughput in ( fT6l ) is plotted against fixed feedback 
delay, which accounts for both CSI quantization and feedback delay. For simplicity, the variable delay 
due to the feedback protocol is omitted but addressed in the sequel. In addition, several capacity metrics 
discussed in Section IVI-Bl are illustrated in Fig. [5] For this figure, the number of antennas in an array is 
L = 4, the SNR is 10 dB, the normalized Doppler shift is //jT = 10~^, and the quantizer codebook size 
is = 128. The quantization loss of 0.2 b/s/Hz is small as also observed in [4], [5]. This loss can be 
controlled by adjusting the codebook size N . The quantization loss is smaller than the capacity loss due 
to feedback delay as it exceeds about 30 samples. In other words, large feedback delay is a dominant 
source of CSI inaccuracy. Finally, the feedback throughput gain decreases exponentially with feedback 
delay. Thus, it is important to constrain the delay for retaining the capacity gain of transmit-beamforming. 

Fig. [6] compares the normalized feedback throughput gain IS.R{D,K) and its approximation in (l26l) 
for different combinations of the fixed feedback delay D and the feedback interval K both in samples. 
Fig.[6ta) and Fig. |6l;b) correspond to the normalized Doppler shift of f£)T = 10~^ and 10^^, respectively. 
The CSI codebook size is = 128 and the SNR is 10 dB. The exact and the approximated values of 
/S.R{D, K) are plotted by using solid and dashed lines, respectively. As observed from Fig. [6l AR{D, K) 
is accurately approximated by ( |26l ). This allows the use of ( [261 ) for computing allowable feedback delay. 
Moreover, AR{D, K) is observed to decrease with D closely following the exponential rate as specified 
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(a) Normalized Doppler shift /dT = 10"" (b) Normalized Doppler shift foT = 10"^ 

Fig. 6. The normalized feedback throughput gain and its approximation versus different combinations of the fixed and protocol 

delays; the normalized Doppler shift is (a) joT = 10"* and (b) joT = 10"^, and the CSI codebook size is Af = 128. 

in Proposition ID This rate grows with increasing normaUzed Doppler shift. This observation agrees with 
the intuition that the effect of CSI feedback delay on throughput is more significant for fast fading and 
vise versa. 

D. Feedback Compression 

Fig. |7] plots the feedback compression ratio Ai? j in (|28] ) for different combinations of normalized 
Doppler shift //jT and the codebook size N . Fig. Ha) and (b) correspond to the normalized CSI feedback 
bit rate of Rf/Rg = 3 and 5, respectively. The threshold for truncating low-probability channel state 
transitions is e = 0.1 (cf. Section IVI-AI ). As observed from Fig. |7J the feedback compression ratio Ai?j is 
independent of normalized Doppler shift fj:,T for the following reason. Roughly speaking, a given value 
of Rf/Rs fixes the degree of temporal correlation in feedback CSI, and thereby results in a constant 
Ai?j independent of JdT. By comparing Fig. Ha) and Fig. Hb), the feedback compression ratio Ai?j 
increases as the normalized CSI feedback bit rate Rf/Rs rises. For example, for the codebook size of 
N = 128, Ai?j is about 15% for Rf/Rs = 3, and rises to about 22% for Rf/Rs = 5. The reason for the 
above observation is that a larger value for Rf/Rs strengthens temporal correlation in feedback CSI and 
thereby increases ARf. Therefore, there exists a trade-off between the feedback compression efficiency 
and the CSI feedback bit rate. 
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(a) Normalized CSI feedback bit rate Rf/Rs = 3 (b) Normalized CSI feedback bit rate Rf/Rs = 5 

Fig. 7. Feedback compression ratio versus different combinations of Doppler shift and the CSI codebook size; the normalized 

feedback bit rate is (a) Rf/Rs = 3 and (b) Rf/Rs = 5, the threshold for truncating low-probability channel state transitions 

is e = 0.1. 

E. Design Example 

This section presents a design example that demonstrates the joint application of the results in this 
paper for designing limited feedback beamforming. This example targets a broadband MIMO downlink 
system employing transmit-beamforming and orthogonal frequency division multiple access (OFDMA) 
[57], which have been included in both the IEEE 802.16 (WiMax) [17] and the 3GPP-LTE [16] standards. 
As illustrated in Fig. [H downlink OFDMA partitions a piece of radio spectrum into frequency slots, called 
subchannels, and assigns them to different users for receiving data from a base station. For simplicity, 
scheduling is omitted in this example and users assigned to subchannels are assumed to be fixed regardless 
of their CSI. Modulated using orthogonal frequency division multiplexing (OFDM), subchannels are 
fully decoupled and each contains a number of orthogonal finer frequency channels, called subcarriers. 
Based on limited feedback, transmit beamforming is performed on each subcarrier to increase its data 
rate. Assume that the narrow-band MIMO channels of different subcarriers in the same subchannel are 
approximately identical. Thus, only one CSI feedback link is required for each sub-channel. CSI is 
estimated at subscriber units by using pilot tones located at the centers of subchannels. 

The design specifications are summarized as follows. The spectrum is 10 MHz at the 2.5 GHz frequency 
band to be used for WiMax implementation. As shown in Fig. [8j this spectrum consists of 8 subchannels, 
each has 1 MHz for data transmission. The remaining 2 MHz bandwidth accounts for pilot and guard 
subcarriers. At each scheduled subscriber unit, the estimated CSI is sampled at 1 MHz equal to the data 
symbol rate. Therefore, the sampling interval is T = 1 microsecond. The size of the codebook used 
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Fig. 8. Spectrum of an OFDMA system at the 2.5GHz band. The lOMHz bandwidth is divided into 8 subchannels of 1 MHz. 
Each subchannel is assigned to one user. 

TABLE I 

Design example of a limited feedback beamforming system 





Specifications 






Bandwidth: 


10 MHz 


Subchannel: 


1 MHz 


Antenna array: 


iVt = iV. = 4 


Carrier frequency: 


/c = 2.5 GHz 


Symbol rate: 


= 1 MHz 


CSI codebook size: 


iV = 128 


Normalized feedback throughput gain: 


{50, 60, 70}% Feedback 


overflow probability: 


Po < 6 X 10"^ 


Mobility: 


y < 45 km/h 







Computed Parameters 





(w/o compression) 


(w/ compression) 




Feedback throughput gain (b/s/Hz) 


0.60 


0.72 


0.84 


0.54 


0.65 


0.76 


(Normalized) 


50% 


60% 


70% 


45% 


54% 


63% 


Throughput gain/subchannel (Mb/s) 


0.6 


0.72 


0.84 


0.54 


0.65 


0.76 


Total throughput gain (Mb/s) 


4.8 


5.76 


6.72 


4.32 


5.20 


6.08 


CSI feedback rate/subchannel (kbits/s) 


21 


21 


21 


17.9 


17.9 


17.9 


Sum feedback rate (kbits/s) 


168 


168 


168 


143.2 


143.2 


143.2 


(Normalized) 


3 


3 


3 


2.6 


2.6 


2.6 


Maximum fixed delay (ms) 


0.99 


0.69 


0.43 


0.99 


0.69 


0.43 



for quantizing feedback CSI is = 128. The normalized feedback throughput gain in (1261 ) is required 
to be {50,60,70}%. The mobility is up to 45 km/h. The feedback outage probability is constrained as 
Po < 6 X 10^^. These design requirements are also summarized in Table H 

Based on above design specifications, the system parameters are computed as follows. Consider CSI 
feedback for an arbitrary subchannel. The Doppler shift for the maximum speed of 45 km/h is obtained 
as //) = 104 Hz [2], and normalized by the CSI sampling rate to be //jT = 10~^. Based on the this 
result and the feedback outage constraint, the corresponding normalized feedback bit rate is obtained 
from Fig. |4] to be = 3. By using Fig. [3] and for the symbol rate of 1 MHz and the codebook 

size of = 128, the CSI source bit rate is obtained as i?s = 3 kbit/s and hence the CSI feedback bit 
rate is i?j = 21 Kbit/s. It follows that the feedback interval is KT = 0.33 millisecond. From Fig. |7J 
the feedback bit rate can be compressed by about 14% for = 128 and /dT = 10^^. Therefore, the 
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compressed CSI feedback bit rate is Rj = 17.9 kbit/s. By combining all subchannels, the sum feedback 
rate is 168 Kbit/s without compression and 143.2 Kbit/s with compression. 

Next, the maximum fixed delay for each CSI feedback link is computed as follows. Corresponding to 
the codebook size of = 128 and the normalized Doppler shift of JdT = 10~^, the eigenvalue defined in 
Proposition m is computed to be \/A = 0.9994 by using the stochastic matrix of the channel state Markov 
chain. By substituting this result into ( [26l ). the maximum fixed feedback delay is obtained as Dmax = 
{992, 688, 431} samples, corresponding to the normalized feedback throughput gain of {50, 60, 70}%, 
respectively. Equivalently, -Dmax = {0.99, 0.69, 0.43} millisecond since the symbol rate is 1 MHz. Last, 
the maximum feedback throughput gain for delay-free CSI feedback is computed to be 1.2 bit/s/Hz by 
using ( fTSl ). Therefore, the required normalized feedback throughput gains of {50, 60, 70}% are converted 
into the actual values of {0.6, 0.72, 0.84} bit/s/Hz. This corresponds to the increase in throughput equal to 
{0.6,0.72,0.84} Mbit/s per subchannel, and {4.8,5.76,6.72} Mbit/s over the whole 10 MHz spectrum. 
The system parameters computed above are summarized in Table Jl 

A few remarks are given on the designed parameters as given in Table IH First, the required CSI 
feedback bit rates in Table U are small due to relatively low mobility and large coherence bandwidth (> 1 
MHz). Second, increasing the feedback throughput gain decreases the maximum allowable fixed delay. 
In other words, higher throughput requires faster signal processing and shorter transmission distance. 
Finally, the CSI feedback reduction by compression is moderate since the normalized CSI feedback bit 
rate and hence the correlation in feedback CSI are small. 

VIII. Conclusion 

We have designed limited feedback beamforming systems over temporally-correlated channels. The 
quantized CSI has been modeled as a first-order Markov chain to provide an analytical tool. This model 
has been validated by simulation. Based on this model, the CSI source bit rate has been derived as a 
function of the Markov chain probabilities. Moreover, adopting a periodic feedback protocol, the CSI 
feedback bit rate supported by the feedback channel has been obtained under a feedback outage constraint. 
Next, an upper bound on the feedback throughput gain has been derived as a function of the CSI feedback 
delay and feedback interval. Last, a feedback compression algorithm has been proposed for reducing the 
CSI feedback bit rate by exploiting temporal correlation in feedback CSI. Simulation results have been 
generated based on spatially i.i.d. Rayleigh fading and the Clarke's fading model. From these results, 
the CSI source bit rate has been observed to increase linearly with Doppler shift. The ratio between CSI 
feedback and source bit rates has been found to be insensitive to a change on Doppler shift. The theoretical 
upper-bound on the feedback throughput gain has been shown to be tight, confirming the exponential decay 
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rate of the feedback throughput gain with feedback delay. Moreover, the proposed feedback compression 
algorithm has been observed to achieve a compression ratio more than 20%. Finally, an design example 
for a limited feedback beamforming system has been presented, which demonstrates the usefulness of 
this work for practical applications. 

This paper opens several issues for future investigation. First, the results of this work focusing on 
limited feedback beamforming can be readily extended to other types of limited feedback systems such as 
precoding, multiuser MIMO, and adaptive modulation and power control. Specifically, the same approach 
based on Markov chain theory as used in this paper can be applied by properly defining Markov chain 
states based the type of limited feedback system. Second, channel estimation errors omitted in this paper 
should be considered in future work. In particular, the training data for CSI estimation should be jointly 
designed with feedback CSI for optimizing system performance and minimizing total overhead. Third, 
a rigorous study of channel prediction for coping with feedback delay is useful for designing practical 
limited feedback systems. Fourth, it is important to investigate the required order of the channel state 
Markov chain for accurately modeling CSI for different kinds of channel distributions such as Rayleigh 
and Rician and temporal correlation. Finally, the time variation of Doppler shift typical in the practice 
should be accounted for in future work. In particular, the time-varying Doppler shift changes system 
parameters including CSI feedback bit and allowable feedback delay. 
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Appendix 

A. Proof of Proposition |2] 

Consider the variation of the channel state in the nth feedback interval (n — l)i^ + l < n < fiK and let 
Zfi denote the number of channel-state transitions in this interval. This does not compromise the generality 
because the channel state Markov chain is stationary. Let denotes the initial channel state for the 

nth feedback interval, thus li^n-\)K = ^o- Conditioned on this initial state, the probability for each sample 
path of the channel state in the nth feedback interval, denoted as ^ ^2 • • • — ^ ^A"> can be expressed in 
terms of channel state transition probabilities as Pr(^i ^ ^2 • • • ^ ^ft: Ko) = PtKtK-\Pe.K-ie.K-i " " " ^t-iUr 
It follows that the conditional probability of no transition is 

Pr(Z^, = I = 4) = [P""] • (35) 
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The conditional probability of single channel-state transition is obtained by grouping all sample paths 
containing only one transition. Thus, 



K N 



Fr{Zn = 1 I I^n-DK =h)=Y.Yl [P'''%mPmio[P''-Xio- (36) 



k=l m=l 



From 051) and ^ 

Vt{Z > 1) = 1 - ^ [Y>T{Zn = I li^k-m = ^) + = 1 I hk-i)K = t)] vr£ 



N 



N K N N 



e=l k=l i=l m=l 



N N N TDK ID ^ / TD \k 

= l-^Pu^i-1^1^ p l^[p^J 

1=1 e=lm=l k=l ^ 

= l-^Pee-i-1^1^ p _ p • 

e=i e=i m=i " 

Thus, by combining the above equation and the feedback constraint, the desired result follows. 



B. Proof of Corollary \J} 

The left-hand-size of ([TT]) attains its maximum value of one at ii' = 1. Thus, the result for 5 ^ 
follows. 

Next, consider the case of — > 1 . An upper bound on K is obtained by considering the following 
upper-bound on the left-hand-size of (fTTl) 

idK , (Pmm- Pu)Pme'^e idK , [max(Pmm, P^)^ - min(Pmm, P«)'^]Pm£7r£ 
^^^'^^ P " l^Pee^i+X. ZT^i 

, , mm J « -1 K mm ^ U\ 

1=1 m^i t=l mj^t 



maXu Paa)^ P-n 
\Pmm Pll\ 



1=1 m^l 



,„ ,Pmm Pu 
m=f=t 



Therefore, the following constraint is the relaxation of that in ([TTI ) 

(max P^af I 1 + X ^r^^^ I > 1 - ^- (37) 

\ ^,\Pmm-Pu\J 

From the above constraint, an upper bound on K, denoted as K^, is obtained as 

log(l - 5) - log (l + Zm^i 

K+ = p , ' • (38) 

log (maXaPaa) 
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To obtain an lower bound for K, the left-hand-size of (fTTI) is bounded from below as 

f;p,fvr,+ j:^%l^l^ > X:^.f-,>(minPj^. (39) 

Thus, the constraint (miiia Paa)^ > 1 — J is more stringent than that in ([TTI l. Based on this constraint, 
a lower-bound on K, denoted as K^, is obtained as 

lofffl - 6) 

^~ = 1 f ■ J V (40) 
Combining (l40l ) and (1381 ) gives the desired result. 

C. Proof of Lemma^ 
From (fTSl) 

Af N 

R = Y^Y.^ [log2(l +P||H„Qj(H„„^)||2) I H„ G V,,H,_d G V^] Pr(H„ G V,,H„_^ G V^) 
£=1 m=l 

Af AT 

II ) I H„ G V^,H„_d G Vm] Pi'(H„ G V^,H„_rf G Vm) 

^=1 m=l 

N N 

= YY ^em Pr(/n = ^, = rn) 

i=l m=l 

where Rim is defined in (ITtI ). Vm is the mth member of the codebook J^, and is the Voronoi cell. 

D. Proof of Corollary \2\ 

By substituting D = into (f76l ) 

AT 

P(0) = [log2(l + P||H„Q/(H„)f ) I H„ G V,] vr,. 

Given that U^V^ = C^"^^^', the result in ([T8] ) follows from the above equation. Next, by substituting (|4]) 
into ([T6l ). the first equality in ( fT9l ) is obtained. The second inequality in ([T9l ) follows from (|7]l, ([TT]) and 
the fact that U^Ve = C^"^^'. 

£. Proof of Proposition |?] 



From (1201 ). the feedback throughput gain is upper-bounded as 

N N 



^=1 m=l 
N / s, N 



iV / \ iV 

< maxP,^ j;|[P^],^-^,| (41) 

m=l ^ ^ f=l 



The desired result is obtained from (|4TI) and Lemma |3] 
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